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The status of lattice calculations of some phenomenology of heavy quarks is presented. Emphasis is on progress 
made in calculating those quantities relevant to estimating parameters of the quark mixing matrix, namely leptonic 
decay constants, the bag parameter of neutral B mixing, and semileptonic form factors. New results from studies 
of quarkonia are highlighted. 



1. INTRODUCTION 

Lattice QCD offers the best hope for estimating 
the non-perturbative QCD effects in weak decays 
of hadrons, thereby allowing the determination 
of elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix from experimental data. Much 
of the phenomenological interest is in those weak 
decays which contain at least one heavy quark, 
but it is this case where technical issues are often 
most difficult and conceptual points most subtle. 

There are a variety of approaches taken to sim- 
ulate a heavy quark on today's lattices. One is to 
use the improved actions used for light-quark cal- 
culations, taking care to keep am not too large, 
and then to rely on an extrapolation (often guided 
by the result of a static calculation) in the heavy 
quark mass from the charm region to the bot- 
tom. A philosophically-different approach is that 
of non-relativistic QCD (NRQCD), namely, to 
integrate out the heavy quark producing an ef- 
fective action which is then discretized; then am 
must be kept large. A third approach, developed 
at Fermilab, seeks to produce an action which 
can be used for any value of am, with systematic 
errors that can be identified and controlled. An 
approximation to the full program is often used 
wherein standard relativistic actions are reinter- 
preted non-relativistically. Much of the debate fo- 
cuses on to what extent systematic errors are con- 
trolled in this case. J. Sloan gave a nice compare- 
and-contrast summary of NRQCD and Fermilab- 
type actions at Lattice '94 

I review developments made and results an- 
nounced this past year. Sec. 2 summarizes re- 



cent calculations of the leptonic decay constants, 
notably fg, where there has been considerable 
progress. In Sec. 3, I present new results, and 
new analysis of old data, for the Bb parameter 
of neutral B mixing. Sec. 4 outlines the suc- 
cesses and limitations of calculations of semilep- 
tonic form factors for B and D decays; this year 
has seen the initiation of some new large-scale 
projects. Sec. 5 highlights new developments 
and results of calculations, presented at this con- 
ference, of the phenomenology of quarkonia. I 
conclude with a qualitative overview of where we 
stand, and in what direction future work may be 
focused. 



2. LEPTONIC DECAY CONSTANTS 

It has been a banner year for the calculation 
of /s, the leptonic decay constant of the B- 
meson. Several large groups, having invested 
many person- years of effort, have presented fi- 
nal results this year; many of these calculations 
have done a comprehensive job of estimating sys- 
tematic errors. There appears to be a consensus 
among many of the groups (See Figure |l}) This 
is encouraging, given the diversity of actions and 
analyses. See reference [g| for a discussion of dis- 
cretization errors. 

The MILC collaboration has published their fi- 
nal results this year ||. They use Wilson quarks 
for both the light and heavy quarks, but use 
some aspects of the "Fermilab interpretation" [Q , 
namely the use of non-relativistic normalization 
and kinetic, rather than pole mass, for the heavy 
quark. MILC has also found it useful to include 
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Figure 1 . A summary of recent lattice results for 
Jb- Data and labels are taken from Table 1. The 
error bars with the longer feet are statistical. Sys- 
tematic errors have been combined in quadrature. 
Also shown is a world average over the subset of 
results which have final and complete estimates 
of systematic errors. 



static- Wilson results to aid in the extrapolation 
in the heavy quark mass. With Wilson heavy 
quarks in the Fermilab interpretation, there is 
a 0(1/ 'M) error (physical, not a lattice artifact) 
which cannot be extrapolated away. MILC, fol- 
lowing JLQCD ||, estimates this effect to be 
1=3 2% from a tree-level approximation. MILC has 
been very careful in estimating systematic errors; 
changes since their preliminary Lattice '97 results 
include updates of their chiral extrapolation, the 
method of combining errors within the quenched 
approximation, and some new unquenched data 
to estimate quenching errors. 

In T. Onogi's Lattice '97 review 0, he reported 
that although the MILC and JLQCD raw data 
agreed for common parameters, there were some 



residual differences from the way the scale was set 
from f n . These discrepancies were discussed after 
last year's symposium by members of each group, 
and the discrepancies have now been resolved |(| . 

The Fermilab (FNAL) collaboration have pub- 
lished their final results [||, updated since Lat- 
tice '97. They used the SW (clover action) for 
both heavy and light quarks, but interpreted 
non-relativistically; this includes the use of cor- 
rectly normalized fields, an additional three- 
dimensional rotation for the heavy quark in the 
current, and the use of the kinetic mass to set the 
quark mass. They include the full 0(a) correc- 
tion, but not the complete O(aa) correction to 
the current. The advantage of their program is 
that they can "sit on" the b mass. They state 
that their systematic errors are under control, 
are smaller than their systematic errors, and are 
smaller than the quenching error which they es- 
timate to be about 10%. 

The GLOK Collaboration has published their 
final results updated since Lattice '97 (see 
also J§|]). Their's is a quenched, (3 = 6.0, calcula- 
tion which uses clover light (tree-level tadpole im- 
proved) and NRQCD heavy quarks. The Hamil- 
tonian includes all 0(1/M 2 ) corrections, the lead- 
ing 0(1/M 3 ) correction, and discretization cor- 
rections at the same level. One important in- 
gredient in the calculation was their axial cur- 
rent renormalization This was the first use 
of mass-dependent matching factors and O(aa) 
current correction. They found mixing between 
lattice currents which does not vanish as M — > oo 
because of an 0(aa) lattice artifact term, which 
can be absorbed into an 0(aaAQco) discretiza- 
tion correction. The latter is anomalously large 
(they found « 12% reduction at the B mass). 
Furthermore, it cannot be included consistently 
if Wilson light quarks are used, that is, if 0(a) 
corrections are included for heavy quarks, they 
must be included for the light quarks as well 
- one is thus left with sizeable scaling viola- 
tions if one uses NRQCD-heavy with Wilson-light 
quarks. GLOK has augmented their (3 = 6.0 cal- 
culation with a [3 = 5.7 run in order to assess 
scaling violations |fTo|| . 

Figure || compares these NRQCD results of Q 
with those of Fermilab W . For the B meson sys- 
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tern, the actions are very similar, 
ought to agree and they do. 



The results 
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Figure 3. 1/M B dependence of $ B from NRQCD 
(for two choices of q*) and Heavy-Clover at (3 = 
6.1. From [Ol. 
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Figure 2. A comparison by Fermilab ||] of their 
results with the NRQCD results of GLOK §. 



JLQCD |Tlj] has presented preliminary results 
at this conference, using actions and analysis 
which follow GLOK's lead. They use quenched 
Wilson gauge actions at (3 = 5.7, 5.9 and 6.1, SW 
light quarks with a clover coefficient calculated at 
one-loop, and NRQCD heavy quarks, including 
both 0(1/M) and 0(1/M 2 ) terms. Their results 
confirm the features seen by GLOK; namely, the 
contribution from operator mixing of the tempo- 
ral component of the axial vector current is larger 
than one would expect || , and this softens the in- 
dependence of Jb, the q* dependence is large at 
one-loop, but is expected to be small at two-loop, 
and the one-loop correction softens the slope of 
Jb versus I/Mb. Figure || shows their compar- 
ison of these preliminary clover-light -NRQCD- 
heavy data with their clover-light-clover-heavy 
data presented last year ||. They are consistent. 

At this conference, Lin and Lellouch jlj] pre- 
sented preliminary high-statistics results of lep- 
tonic decay constants using a tadpole-improved 
tree-level SW quark action for both heavy and 



light quarks (with KLM normalization for the 
quark fields). They ran at j3 = 6.0 on a 16 3 x 48 
lattice and at /? = 6.2 on 24 3 x 48. They do not 
attempt a continuum extrapolation, so I, as they, 
will take their (3 — 6.2 data as their best. 

The ALPHA collaboration |b|] has proposed a 
method for systematic improvement of the Wil- 
son action and bilinear quark operators to re- 
move all 0(a) discretization errors. The program 
evaluates the clover coefficient and renormaliza- 
tion constants (of bilinear quark operators) non- 
perturbatively. This year, two new projects have 
calculated heavy-light leptonic decay constants 
using this approach. APE Q and UKQCD (T|| 
ran at = 6.2 on 24 3 x 64 and 24 3 x 48 lattices, 
respectively. Both simulate with the heavy quark 
near charm to keep 0(am) errors small; thus they 
obtain jn by interpolation. However, APE finds 
that the extrapolation up to the bottom region is 
quite sensitive to the functional form of the fit; 
this inflates their systematic error. In these pro- 
ceedings, UKQCD quotes only statistical errors. 

Tables 1 and 2 list and Figure [l] displays se- 
lected results for heavy-light decay constants and 
ratios. In computing a world average (Table 3), 
always a notoriously difficult thing to attempt, I 
include only those final, not preliminary, large- 
scale results where modern methods have been 
used and a comprehensive investigation of sys- 
tematic errors has been done, which includes ei- 
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177(19)(S) 
L 162(?)(?|) 
12] 161(16) 
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190(12) 


1.16(4) 

1.18(3)(|) 
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147(ll)(f 2 )(9)(6) 
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173(4) (12) 
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171(10)(f)(f) 

175(8)(I )(11)(7)( 7 ) 

185(P)(9) 

199(3)(14) 
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i-n(M)(l)(l) 

1.20(4)($) 
1.13(1) 
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PCW93 |7J 
UKQCD93 |[ 
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180(50) 
1 160(g)(f|) 

187(10) (34) (15) 


194(l)(f) 
207(9) (34) (22) 


1.09(2)(5) 

1.22(f) 

1.11(6) 



Table 1 



Results for fs, fs s and JbJ Jb obtained in the quenched approximation. The errors are statistical, 
then various systematic. Results are divided into three sections. In the top section are the most recent 
preliminary results; thus, despite their quality they are not included in world averages. In the middle 
section are recent published and finalized results where all systematic errors (except quenching) have 
been carefully assessed; these are folded into a world average. In the last section are selected older results 
from pre-modern methods. 



ther a continuum extrapolation, or an estimate of 
scaling errors. Also tabulated and graphed are a 
sampling of other results. Emphasis is on recent 
results, especially those preliminary results pre- 
sented at this conference. See reference 0] for a 
more comprehensive listing of older results. 

Also reported at this conference were two inves- 
tigations into the spectrum of heavy-light mesons. 
Ali Khan [^2j reported on results from an exten- 
sive calculation of B, B c T, and 6-baryon spec- 
tra; the study included radially and orbitally ex- 
cited mesons. Lewis and Woloshyn 23 studied 



the relative effect of 1/M, 1/M 2 and 1/M 6 terms 
from the NRQCD heavy quark. The expansion 
is well behaved for bottom; effects are larger, as 
expected, for charm. 

In summary, leptonic decay constant calcu- 
lations are quite mature; several large groups 
and years of effort have resulted in results with 
small statistical errors, and systematic errors 
which, aside from quenching, are largely iden- 
tified and controlled. Pleasingly, there appears 
to be a consensus on values among those groups 
that have recently published their final results. 
(This consensus needs to be confirmed by the 
new projects which use non-perturbative esti- 
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Table 3 

World averages for quenched f B , /s a , Id, Id s , 
and ratios. The "unquenched" results are best 
guesses obtained by shifting the quenched val- 
ues by the rough estimates of quenching errors 
from MILC Q (the last of their errors in Tables 1 
and 2). Estimates of unquenched errors assume 
that the MILC corrections could be off by a factor 
of 2. 
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Table 2 

Same as for Table 1 but for /d, fr>„, !dJ !d- 



mates of coefficients and extrapolate up from 
lighter masses.) The world average of fs in the 
quenched approximation, obtained from from var- 
ious actions and groups, has stabilized in last 2-3 
years (world averages had been reported [p|,^4|, 
llJIHQ in the range 160-175 MeV with errors 
of about 25-35 MeV), and is significantly lower 
than several years ago (when world averages were 
200(40) MeV [§7| or above). Looking back @, 
one sees that older results were higher as sev- 
eral effects conspired: old static values which 
guided extrapolations in 1/Mb were misleadingly 
high prior to modern variational techniques which 
remove excited-state contaminations, and both 
1/Mb corrections to the static limit and dis- 
cretization errors are larger than were expected. 
Quenching errors currently are difficult to esti- 
mate. Unquenching may increase fs by 10-15% 
(MILC), maybe more p8fl, or maybe less p9| . 
Lattice estimates of fs continue to be very useful 
phcnomenologically. In the future, we'll see con- 
tinued work with improved actions and renormal- 
ization, with the emphasis shifting towards mak- 
ing better estimates of quenching errors. 

3. NEUTRAL B MIXING 

The CKM matrix element |Vtd| can be ex- 
tracted from the mass difference induced in neu- 
tral B mixing 

(Am)| B _ 5 = f 2 B B B \V td \ 2 F(m t ,M w ) 



where F is a known function and the Bb param- 
eter has been used to parameterize the matrix 
element (B°\O L \B ) of the AB = 2 four-quark 
operator O l = ^(1 - 7s)<7 &7m(1 _ 7s)? in terms 
of its approximation under the factorization hy- 
pothesis ("vacuum saturation" or VSA) 

Bb = (B°\O L \B ) = (B°\O L \B ) 
(B°\O l \B°) vsa 8/3/|Af| 

Since one needs a non-perturbative evaluation 
of the matrix element and thus the combina- 
tion f B Bs, why calculate Jb and Bb separately? 
Reasons include (1) a quite precise value can 
be obtained, with an optimal choice of smear- 
ing function, from relatively few configurations, 
because B B can be extracted from a ratio of 
three- to two-point functions which are strongly 
correlated, (2) perturbative corrections tend to 
be stabilized because of cancelations in numera- 
tor and denominator, and (3) it seems as though 
VSA may be a surprisingly good approximation 
(to k 10-15%) for the B B parameter; this is 
an important qualitative statement, of use to 
model builders, which should not be obscured by 
poor-statistics attempts to calculate the product 
b bSb- 

The Bb parameter was reviewed by J. 
Flynn ||| at Lattice '96 and by A. Soni [|o| at 
Lattice '95. Results from UKQCD (Lin and Lcl- 
louch) and from the Hiroshima group fl3l| |, 
were reported at this conference. Other recent 



G 



results include those from Di Picrro and Sachra- 
jda from Gimenez and Reyes P^j , and from 
Bernard, Blum and Soni j34J]. 

This year, Lin and Lellouch have calcu- 
lated the Bb parameter as part of larger effort to 
calculate B parameters and decay constants for 
light and heavy mesons. They use the SW ac- 
tion, with the clover coefficient tadpole improved 
at tree-level, at both (3 = 6.0 and 6.2 to check 
scaling. They see little a-dependence. In fact, 
the world's collection of results from conventional 
(non-static) methods (Figure ||) (the top set of 
results) shows consistency among groups, no a- 
dependence and agreement between Wilson and 
SW calculations. 
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Figure 4. Heavy quark mass dependence of Bb 3 - 
From UKQCD (Lellouch/Lin) @. 

Figure ^ shows the heavy quark mass depen- 
dence of Bb s - There is a consensus Jl^Q] that 
the extrapolated-static value is larger than that 
at the physical mass. 

The extrapolated-static results and those ob- 
tained directly in the static approximation should 
agree. There are disagreements among the static 
calculations, some of which do not agree with 
extrapolated-static results (Figure ||). As with 
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Figure 5. World results for the BB(rrib(mb)) pa- 
rameter. For a more complete listing, see ta- 
bles in |35| and pl|,p6| ]. "Conventional meth- 
ods" by UKQCD |L2| , JLQCD (36), ELC @, 
and Bernard & Soni @ labeled "Clo-Clo" or 
"Wil-Wil" use SW (clover) or Wilson, respec- 
tively, for both heavy and light quarks; the world 
average is taken from among this set. "Ext-Clo" 
and "Ext-Wil" are the extrapolations of these 
data to the static limit for the heavy quark. 
To these should be compared results of using 
directly the static approximation for the heavy 
quark (from Kentucky [ff], UKQCD || and 
Gimenez & Martinelli [3Jj]) which are labeled 
"Sta-Clo" for a clover, or "Sta-Wil" for a Wilson, 
light quark. Gimenez and Reyes p3| have reana- 
lyzed the renormalization of the old data of refer- 
ence |H and of reference jQ (labeled UKQCD*). 
Di Pierro and Sachrajda ||32(| have made a similar 
reanalysis and reach the same conclusions. 
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decay constant calculations, one must fight the 
very poor signal-to-noise ratio inherent with using 
static quarks. A large number of configurations 
has to be generated if one uses canonical smearing 
techniques [|9| . The Kentucky group obtains |55| 
the same statistical error with an order of magni- 
tude fewer configurations by using optimal smear- 
ing techniques that they developed Q ■ It is no- 
table that all groups' raw lattice data are con- 
sistent; furthermore, the VSA holds surprisingly 
well for all operators (with appropriate normal- 
ization). 

As a short digression, it is interesting to note 
that factorization seems to hold well in related 
static calculations. Di Pierro and Sachrajda |32j 
need to evaluate matrix elements like (B\ &7 M (1 — 
75)? <77/j(1 — 75)^ \B) for their studies of inclu- 
sive decays of heavy hadrons. The operators are 
similar to, but different than, those for the Bb 
parameter. For the "figure-8" contractions, they 
find that factorization holds surprisingly well — 
to w 5%. 

For the static Bb parameter calculations, dif- 
ferences in renormalized quantities arise from 
choices made in lattice-to-continuum matching: 
whether products and quotients are expanded in 
a or merely multiplied and divided makes a un- 
comfortably large difference. 

Although all organizations of perturbation the- 
ory at one-loop are theoretically equal, some 
arc more equal than others! The Kentucky 
group |35| has been advocating an organization 
of lattice perturbation theory which uses the 
Lepage-Mackenzie Q choices for a s and tadpole 
improvement and which, by allowing explicit can- 
celations, is insensitive to the wave-function nor- 
malization and which does not mask the agree- 
ment of the raw data with VSA. (This also re- 
duces the statistical errors.) Their results agree 
with VSA and with the extrapolated static re- 
sults and can thus be used in conjunction with 
conventional data to interpolate to the B-meson 
mass, as has been done with much success for the 
decay constant calculations. The original clover- 
light-static-heavy results |38|]39f| which disagreed 
with the Kentucky results and with the world's 
extrapolated-static results, have recently been re- 
analyzed [|32|]33j| with an organization of perturba- 
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2151S 




1.14±ij 


1.1415 



Table 4 

World averages for quenched and "unquenched" 
quantities. 



tion theory which includes tadpole improvement 
(and corrects errors in some renormalization con- 
stants); the discrepancies have been reduced but 
not eliminated. 

The Hiroshima group |3lJ] has begun the first 
calculation of the Bb parameter using NRQCD 
for the heavy quarks (with 16 3 x 48 quenched Wil- 
son j3 = 5.9 gauge-field configurations, and SW 
light quarks with a tree-level tadpole-improved 
clover coefficient). They can see effective-mass 
plateaus despite using local-local correlation func- 
tions. They find large dependence on the heavy 
quark mass. For their preliminary calculation 
they rely on the static renormalization constants. 
It is not clear that the large mass-dependence will 
be mitigated upon using the correct heavy-mass 
dependent renormalization constants, as was the 
case for the decay constant . 

Relying, then, only on conventional results, I 
quote a world average for Bb(ttib) in Table 4. 
This value has been historically stable; it agrees 
with Soni's Lattice '95 estimate 1 30 . Indeed, 
even prehistoric estimates M2| are good because 



any errors in normalization (e.g. non-KLM fac- 
tors) cancel in the ratio. All groups agree that 
the light quark mass dependence of Bb is small. 

Bernard, Blum and Soni (BBS) ]34| have rec- 
ommended that the SU (3) flavor dependence be 
investigated by calculating directly a ratio of ma- 

trix elements giving » f f s f f 3 „ 3 One would ex- 

pect that many systematic errors would cancel. 
Lin and Lcllouch filial have since calculated this 
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quantity using clover quarks (BBS used Wilson). 
Both groups find that the direct method is com- 
patible with, but unfortunately has larger statisti- 
cal errors (by a factor of 2 or 3) than, the indirect 
method which combines separate calculations of 
Bb and fs, for both strange and light quarks. 

4. SEMILEPTONIC DECAYS 

The differential decay rate for P(Qq) — > 
P'(q'q)lui is 



dT 
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\V q , Q \ 2 \f + {q 2 )f 



where f+(q ' ) is a vector form factor, obtained 
from the matrix element 



(PV)|vgp(p)> 



TOp — TO , 



= [p + p - q- 



nvh — m. 
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Similarly, the differential decay rate for 
P(Qq) — ► P*{q'q)lvi is expressed in terms of 
a vector form factor V(q 2 ) and two axial form 
factors Ai(q 2 ) and A 2 (q 2 ). 

Experiments are able to determine an intercept 
and slope for / + 1 ( ) , so that it is conventional to 
quote theoretical predictions at q 2 = 0; further- 
more, the form factor is often parameterized with 
pole dominance form f + (q 2 ) — f + (0)/(l—q 2 /m 2 ) 
to aid interpolation and extrapolation. Unfor- 
tunately, lattice calculations work best at small 
three-momenta, i.e. near zero recoil, = 
(mp — rnp/) 2 , and so for decays such as B — > 7r, 
large extrapolations are necessary. 

Semileptonic form factors were reviewed by T. 
Onogi at Lattice '97 §|, by J. Flynn f|| at Lat- 
tice '96 and by J. Simone |^3| at Lattice '95. For 
a very comprehensive and recent survey of lat- 
tice results, see the review by Flynn and Sachra- 
jda (on which the discussion of the present 
section relies heavily). 

4.1. Semileptonic D — * K,K* Decays 

For semileptonic D — * K, K* decays, the max- 
imum meson recoil momenta are below 1 GeV/c 



so all q 2 can be sampled with small momentum- 
dependent discretization errors; thus one can in- 
terpolate to q 2 — 0. Furthermore, one can sit on 
s and c so the only mass-extrapolation needed is 
for the light quark. Thus, D decays are an excel- 
lent test for the lattice; experimental determina- 
tions of \V CS \ are already good and are likely to 
improve soon. Some disadvantages of the lattice 
calculations are that the charm and strange quark 
masses are non-degenerate, so there is no natu- 
ral normalization condition, and that care need 
be taken to reduce discretization errors for the 
charm mass. 



Experiment 



Clo-Clo_APE95 



Clo-Clo_UKQCD95 



Wil-Wil_LANL95 



Wil-Wil_WUP97 




Figure 
APE § 
pertal 



.5 .6 .7 .8 .9 
f + (0) (D^K) 

Lattice estimates of / + (0) from 
UKQCD H|, LANL |§ and Wup- 
, compared to experiment pq|. 



Figure |^ compares a recent experimental esti- 
mate |Q with a collection of lattice results for 
/+(0), from APE @, UKQCD @, LANL g§] 
and Wuppertal ^?]]. (See |2l[] for a more complete 
summary including older results; those listed here 
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use either the SW action, or Wilson with KLM 
normalization.) The agreement is quite satisfac- 
tory. The agreement for the D — > K* form fac- 
tors, V^(0), .Ai(O), and A 2 (0), is also good, al- 
though somewhat poorer for A\ and A2, which 
depend upon correctly normalizing the lattice ax- 
ial current. 

All of these results are from quenched calcu- 
lations without a continuum extrapolation. A 
new effort this year seeks to remedy this defi- 
ciency. Preliminary results from FNAL were pre- 
sented by J. Simone at this conference |^9|. Their 
quenched calculation uses the SW (clover) ac- 
tion interpreted non-relativistically in the "Fer- 
milab approach" ; thus they can "sit on the charm 
quark". At this early stage in their project, the 
light quark has strange-quark mass. They see a 
gentle a-dependence of the matrix elements which 
bodes well for controlling the continuum extrapo- 
lation. Their goal is to precisely determine decay 
rates for D -> nlv and D -> Klv; the FOCUS 
experiment will soon have a high-statistics deter- 
mination of BR(D -> wlv)/BR(D -> Klv). 

Preliminary results from UKQCD were pre- 
sented by C. Maynard at this conference [f50f . 
Their quenched calculation on 24 x 48 lattices 
((3 — 6.2) uses a O(a)-improved SW with a non- 
perturbative value of clover coefficient, and sim- 
ulates with a heavy quark mass near charm, and 
a light quark mass both near strange and close 
to the chiral limit. / + , f° are simultaneously fit 
(to pole-dominance model) for spatial and tempo- 
ral components of vector current. The statistics 
are quite competitive: 3-5% errors for D — > K 
(from 216 configurations). The results are to be 
compared with older results which used tree-level 
clover-coefficients . 

4.2. Semileptonic B — > 7r,p Decays 

Semileptonic decays of the B are very impor- 
tant for constraining CKM matrix elements. Lat- 
tice calculations are really needed since HQET 
symmetry is not as helpful as for B — > D. Un- 
fortunately, the large bottom mass means that 
simulations must be at the charm mass and then 
extrapolated to bottom, else care must be taken 
to design an action for which discretization er- 
rors are small and controlled. Furthermore, the 



maximum meson recoil momenta are large so all 
q 2 cannot be sampled (contrast D decays) with- 
out large momentum-dependent discretization er- 
rors. Thus large extrapolation to q 2 = are made 
for which one must rely on the pole dominance 
model. 

Although heavy quark symmetry says noth- 
ing about overall normalizations, it does predict 
that hadronic matrix elements for light-meson de- 
cay modes such as P — > ttIv, scale in the heavy 
quark limit, M — > 00, for q 2 near q^ax ( zer0 re- 
coil). Form factor scaling follows: the leading 
mass dependence is f° ~ M^ 1 / 2 , /+ ~ M +l / 2 
for B -> nlv, and V ~ M+ 1 / 2 , A 1 ~ A/" 1 / 2 , 
and A 2 - M +1/2 and for B — > plv. There is a 
consensus (see T. Onogi's review [g] at last year's 
lattice conference) that such scaling is seen in lat- 
tice simulations for q^ax' an< ^ thus can be used to 
extrapolate simulated form factors to the B mass. 




Figure 7. Fit to lattice results for V, Aq, Ai, Ti, 
T2, versus (9/GeV) 2 , for p final state. From [fill . 

To obtain results at smaller q 2 , one has to ex- 
trapolate aided by a pole model. UKQCD (|l) 
has reanalyzed some older data and exploited the 
scaling constraints of the light-cone sum rule to 
guide model-dependent extrapolation to q 2 = 
where the leading mass dependence for all form 
factors is ~ M -3 / 2 . (They point out that the 
analyses of some other groups violate these scal- 
ing relations.) Figure [j] show the results of a si- 
multaneous fit which respects all constraints. One 
can see that the extrapolation to q 2 = is very 
large indeed. 
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Wil-Wil_ELC94 




.0 



.6 



.1 .2 .3 .4 .5 
f + (0) (B^TT) 

Figure 8. Lattice estimates of B — » 7r / + (0) 
from UKQCD @, Wuppertal |47), APE @ and 
ELC 



Figure g compares previous lattice estimates of 
/+ at q 2 = from UKQCD Q, Wuppertal J4T) , 
APE Q and ELC ||. 

Preliminary results from UKQCD were pre- 
sented by C. Maynard at this conference 50| , 
in conjunction with their semileptonic D decay 
calculation (quenched, 24 x 48, (3 — 6.2, 0(a)- 
improved SW with a non-perturbative clover co- 
efficient , heavy quark mass near charm). At last 
year's conference, we saw evidence || of a large 
violation of the soft-pion theorem which predicts 
the relation 



2 ) 

max / 



/b//tT 



This year, UKQCD sees no such violation as 
demonstrated in Figure §. The UKQCD cal- 
culation includes renormalization constants, and 
uses a smaller lattice spacing [/? = 6.2, versus 
f3 = 5.9 for JLQCD || (SW action) and /3 = 5.8 
for Hiroshima [tl (Wilson light quarks, NRQCD 



F vs 1 /M H q max 
□ F=f B /f„ (w) 

o F=fo(qLJ M 



2 — • 



1 ' M 



Figure 9. 
From flscll 



Comparison of /°(g max ) and /b/U 



heavy)]. 

FNAL, as presented by S. Ryan at this con- 
ference, are engaged in a multi-a calculation of 
B semileptonic decays, in conjunction with their 
semileptonic D decay calculation, with the inten- 
tion of taking the chiral and continuum limits. As 
described in the previous subsection, the advan- 
tage of using their action is that they can sit at 
the bottom mass (or any other), and avoid hav- 
ing to extrapolate in the heavy quark mass. They 
verify that the heavy-mass dependence of the 
matrix elements is very gentle, from the charm 
through the bottom regime. Of course, they are 
subject to the same limitations as other groups 
for extrapolations in the light quark mass, and in 

A novel approach |5^] has been proposed to pre- 
dict form factors for exclusive processes such as 
B — > 7r by computing light-cone wave functions 
(not, as before, just their moments) in terms of 
lattice correlation functions. Although no simula- 
tion results have been presented yet, the method 
is intriguing because it can be used at small q 2 
directly. 

In summary, all groups agree on the 1/M scal- 



ing near 



ELC, APE and UKQCD use 
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scaling to extrapolate from charm to bottom; 
JLQCD, Hiroshima and FNAL can sit on the 
bottom. Estimates of the form factor at q 2 = 
have a large extrapolation, which is model depen- 
dent. Continuum extrapolations have not been 
published, but are in the works. A comprehen- 
sive estimate of systematic errors, as for fs, has 
yet to be made. One can hope that future esti- 
mates of systematic errors might follow improve- 
ments in decay constant calculations, with a lag 
of two or three years. It will continue to be dif- 
ficult to simulate at low q 2 , but it may not be 
necessary. Experimental data may soon be avail- 
able at high q 2 ; if so, a comparison with accurate 
lattice results will be sufficient to determine \V u b\ 
with impressive precision. 

4.3. Semileptonic B — » D, D* Decays 

Semileptonic decays B — > D,D* are used to 
extract the \V c b\ CKM matrix element. Heavy 
quark symmetry plays a central role. The intu- 
itive picture is well known. (See reference [f56f , 
for example). Typical momenta exchanged be- 
tween heavy and light quarks are O(Aqcd); 
light degrees of freedom cannot resolve distances 
1/toq <C 1/Aqcd and are blind to the mass and 
spin of heavy quark; thus hadronic systems which 
differ only in the mass and spin of the heavy quark 
have the same light dof configuration. The non- 
perturbative light dof ("brown muck") interacts 
with heavy-quark static color field (which extends 
over long distances) . When an external weak cur- 
rent boosts the heavy quark (rat — > ra c ), the 
brown muck reacts and the form factor suppres- 
sion is universal 

Heavy quark symmetry suggests rewriting the 
matrix elements in terms of form factors which 
are functions of velocity transfer lu = v-v , rather 
than momentum transfer. Thus (P' | P) oc 
(v + v')h + (uj) + a(v — v')h-(uj). Similarly, 
(P*|V^|P) is expressed in terms of hy and 
(P*\A fl \P) in terms of h^, Iia 2 , ^a 3 - Then 

hi(u)) = [oti + f3i(m b ,m c ,Lu) + 0(A/toq)] 

where the HQS symmetry enforces a.+ t v,A u Aa — 1 
and a_ ; A 2 — @- C( w ) IS the famous Isgur-Wise 
function. 

The experimental measurement of the differen- 







Pt \V cb \x 10" 


BSS |6| 
UKQCD_p| 
LANL g§ 


o.9(!)(!) 

0.97(6) 


1.4(2)(4) 4.4(5)(7) 
1.2(I)(f) 3.7(1)(2)(J) 



Table 5 

A collection of lattice results for the slope of the 
Isgur-Wise function. This is then used to con- 
strain extrapolation of the experimental data to 
obtain the CKM matrix element listed. 



tial cross sections yields |T4i,| 2 J-"(ci;) 2 where the 
physical form factor, T(u>), equals £(cj) plus per- 
turbative and power corrections. Experimentally, 
the shape of !F(uj) is not well known, and so it 
is common to parameterize the extrapolation of 
the data back to zero recoil (where a theoreti- 
cal calculation of the normalization .F(l) would 
then determine in terms of the slope at 

zero recoil, — p 2 . A recent world (experimental) 
s2 



average is p" — 0.75(11) |57|. A more precise 
theoretically-determined value would greatly aid 
the extrapolation back to zero recoil. A task for 
lattice calculations is then a determination of p 2 , 
or of the corresponding slope, — p 2 , of the Isgur- 
Wise function to which it is related by calculable 
corrections Q. 

The first lattice calculations of the Isgur-Wise 
function were reviewed by Kenway at Lattice '93. 
(See also the aforementioned reviews by J. Si- 
mone Q and by Flynn & Sachrajda |pl[.) Since 
then, there have been surprisingly few lattice cal- 
culations, some of which are listed in Table 5. 
There also exist calculations p8|-|6"o[ which simu- 
late the heavy quark effective theory directly on 
the lattice, but renormalization is subtle, and dis- 
cretization errors can be sizeable for lowest-order 
actions on modest-sized lattices. 

There has been renewed interest recently, and 
at this conference we had reports of simulations 
from FNAL (presented by S. Hashimoto p3|) and 
UKQCD (D. Richards pi). A. Kronfeld|[64) , P. 
Boyle |65| and J. Sloan |6(| reported on renormal- 
ization calculations. 

The UKQCD calculation is part of a larger cal- 
culation including D — > K, K* and B — > w, p form 
factors as described above. Recall that the clover 
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coefficient is determined non-perturbatively. The 
matching coefficient in the renormalized vector 

current = Z v (l + b v am q )[V^ att H ] can be 

computed for degenerative transitions at zero mo- 
mentum, and then compared with the ALPHA 
collaboration's |l3| non-perturbative determina- 
tion of Zy and by . The striking agreement leads 
UKQCD to believe that they have good under- 
standing of the discretization errors. 

S. Hashimoto (6^] presented preliminary FNAL 
results for B — > D, using the same action as for 
the other semileptonic form factor calculations 
previously described. Recall they can simulate 
at the b mass; currently, the light quark mass is 
around strange. They consider a ratio of matrix 
elements at zero recoil 



\h+(l) 



(D\cj°b\B){B\b-/°c\D) 
{D\cj°c\D){B\b~/%\B) 



for which statistical errors are small, current 
renormalization and systematic errors largely 
cancel. (This is the same ratio as first defined 
by Mandula and Ogilvie |58| and then used by 
Draper and McNeile |fj§ for lattice HQET.) Us- 
ing a corresponding ratio to obtain they then 
construct a linear combination to obtain ^(1), 
i.e. the physical form factor at zero recoil! Previ- 
ously, lattice calculations had predicted only the 
slope to constrain the fit to experimental data 
leaving ^(l^V^I as the free parameter, and then 
relied on sum rule calculations of .F(l). With this 
development, HQET and the lattice alone can de- 
termine \V c b\- 

5. HIGHLIGHTS FROM QUARKONIA 

At this conference we saw several interesting 
calculations of quarkonia. I will not attempt to 
put these results in context by making a compar- 
ison with a complete list of earlier results, due to 
space constraints, but rather will take this oppor- 
tunity to advertise several of these calculations. 

Presently, the usual approaches used in calcu- 
lating properties of quarkonia on the lattice are 
NRQCD (expansion in v 2 ; breaks down for small 
masses and small lattice spacing) and the Fermi- 
lab approach (improve at am 3> 1; coefficients are 
mass-dependent; present simulations with "non- 



relativistic reinterpretation" of Wilson-SW ac- 
tion). 

Problems show up in charmonium hyperfine 
splitting. At this conference, Shakespeare and 
Trottier |57) compared the effect of 0(v 6 ) versus 
0(v 4 ) terms, and the effect of the detailed way in 
which mean field is invoked (mean-link from Lan- 
dau gauge versus from fourth-root-of-plaquette 
tadpole improvement). They see NRQCD break- 
down for oM c ° < 1 (as expected) and see a clear 
preference for using the Landau prescription. It 
should be noted that this is not an indictment 
of NRQCD, nor of mean-field improvement. Cal- 
culations of bottomonium are well under control. 
For charmonium, the effect of higher order terms 
is as expected, except that the prefactors are a 
factor of two or so higher than one would predict 
from dimensional analysis. 

T. Klassen showed us the first results from 
a long-term program which provides an alterna- 
tive to NRQCD or Fermilab: Symanzik improve- 
ment on anisotropic £ = a s /a t S> 1 lattices. His 
quark action shares some formal similarities with 
that of Fermilab. It has on-shell 0(a) classically- 
improvement plus a non-perturbative O(a) im- 
provement which requires a tuning of the Wil- 
son parameter r(m, £). However, the coefficients 
of the electric and magnetic terms in the action 
are mass-independent (contrast Fermilab). For 
his first calculation, he uses mean-field (Landau) 
to estimate these coefficients, and obtains r(m, £) 
from charmonium dispersion. His continuum ex- 
trapolation is well-behaved. He sees spectacular 
agreement with experiment for the charmonium 
P-state splitting. 

R. Horgan J69[ presented preliminary results of 
a long-term project using an anisotropic lattice 
for NRQCD. Their goal is a high statistics calcu- 
lation for bottomonium and for hybrid mesons. 
Early results are encouraging and are in good 
agreement with those from isotropic lattices. 

K. Hornbostel j70| presented a new analysis 
from the NRQCD collaboration of the bottom 
(and charm) mass. They used a tadpole-improved 
NRQCD action including next-to-leading order 
relativistic and discretization corrections, com- 
puted spin-averaged IP-IS splittings to deter- 
mine the scale a, tuned the bare mass so the 
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kinetic mass from the dispersion relation agreed 
with experiment, perturbatively related the bare 
mass to the M S mass (using a Lepage-Mackenzie 
q*), and ran the MS mass to its own scale. They 
obtain 4.28(3)(3) GeV from their = 6.0 data 
(which has the best statistics). 

Bali and Boyle |7l]] have done some nice work in 
which they measure the potential on the lattice, 
and then parameterize it to solve the Schrodinger 
equation on (in)finite lattices. They then vary 
the lattice volume and quark masses to deduce 
from their model qualitative features seen in the 
full simulations; namely, non- negligible finite- 
size effects, a non-trivial dependence on quark 
mass, and splittings strongly dependent on the 
Coulomb coupling. 

I would also like to advertise a calculation of 
the heavy-quark leptonic width using NRQCD by 
Jones and Woloshyn [jf2|, a new estimate of a s 
from the SESAM collaboration ]73|] , and a study 
by Fingberg of the effect of the gluon conden- 
sate on the low-lying quarkonium spectrum. Pen- 
nanen [7j| presented results from a study of four- 
quark systems, while Morningstar [[r6[ presented 
final results for the presence of gluonic excitations 
in the presence of a static quark-antiquark pair. 

6. CONCLUSIONS 

The use of improved actions, techniques, anal- 
ysis and computer power in recent years has been 
most fruitful. Lattice QCD has made very useful 
and competitive predictions of many phenomeno- 
logical quantities crucial for the confrontation of 
theory and experiment. 

It has been a banner year for the calculation 
of heavy-light leptonic decay constants such as 
Jb- There are now reliable estimates of system- 
atic errors for the quenched calculation, and there 
is good agreement among groups using a variety 
of methods. In the near future, we should see 
if the final analysis (including a continuum ex- 
trapolation) from conventionally-interpreted SW 
actions (which use a non-perturbative estimate 
of the clover coefficient) gives results which con- 
cur with the clustering of several results obtained 
from non-relativistic interpretation. If so, atten- 
tion must turn to the remaining formidable is- 



sue of more reliably estimating the effects of the 
quenched approximation. With new actions and 
methods, more attention should be paid to the 
heavy-light spectrum, where current estimates of 
the sensitive hyperfine splitting are too small. 

Estimates of the Bb parameter from neutral B 
mixing seem rather robust; since it is determined 
from a ratio of correlation functions, systematic 
errors and statistical fluctuations largely cancel. 
More modern estimates using Wilson or SW ac- 
tions agree with older results. The vacuum satu- 
ration, or factorization, approximation appears to 
hold at the 15% level. Earlier variance among re- 
sults using the static approximation for the heavy 
quark are better understood, and somewhat re- 
duced. Simulations with an NRQCD heavy quark 
were presented this year, but need a difficult 
renormalization calculation for completeness. 

For semileptonic decays, we need need more 
complete estimates of systematic errors (contin- 
uum limit, variety of actions). It would be nice 
to see more groups calculating these form factors. 
The kinematics of D — > K,K* provide a good 
place to test actions, techniques and analysis. For 
B — > 7T, p, the kinematics demand a large, model- 
dependent, extrapolation in q 2 . (One must also 
beware discretization errors due to the large bot- 
tom mass; the two philosophies used are the same 
as for leptonic decay calculations: either attempt 
to control errors with a non-relativistic interpre- 
tation of the action, or simulate around the charm 
mass, and boldly extrapolate.) It is unlikely that 
the large q 2 extrapolation can be avoided in the 
near future with conventional approaches. It may 
not be necessary. We expect our experimental 
colleagues, using increased luminosities available 
around the turn of the millennium, to be success- 
ful in measuring the form factors at high q 2 . Then 
a direct comparison of theory and lattice results 
should provide a good estimate of \V u b\- 

The power of heavy-quark-symmetry applied to 
B — > D semileptonic decays has reduced the phe- 
nomcnologist's dependence on lattice estimates, 
but the lattice can still make fruitful contribu- 
tions by providing a reliable estimate of the slope 
and intercept of the form factor at zero recoil. 
New calculations using relativistic heavy quarks 
have been presented at this conference; it is wel- 
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come to see renewed interest in these form fac- 
tors. We expect NRQCD estimates to be avail- 
able soon and anticipate that within a year or two 
these, relativistic estimates, and estimates from 
directly simulating HQET, will all be reconciled 
as has been the case for leptonic decay constants. 

Calculations of quarkonia continue to be a good 
forum for testing new actions and techniques 
and for predicting parameters of the standard 
model; we saw updated estimates of a s and rrif, 
at this conference. Bottomonium calculations are 
in good shape; NRQCD calculations in particu- 
lar are quite precise. The convergence of the ex- 
pansion in inverse mass is poor for charmomium, 
however, and so this is a good place to try new im- 
proved actions. We saw that anisotropic actions 
look very promising. 

Calculations of heavy quark physics on the lat- 
tice have reached a new level of maturity in the 
last few years, and have provided the wider com- 
munity with reliable non-perturbative estimates 
of several quantities of phenomenological inter- 
est. But there are still many quantities for which 
it is difficult to be as precise. Meanwhile, new 
approaches and techniques mean that the lattice 
calculations will continue to interesting and im- 
portant for years to come. Onward and upward! 
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